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(54) A thermal design method and apparatus 

(57) A design method for structures is invented, tak- 
ing into ac count the effects of fluid temperature fluctu- 
ations, comprising: numerical calculations of time-aver- 
aged solutions for the motion equation and energy 
equation for fluids; a Process 1 in which the low-cycle 
components of the temperature fluctuations occurring 
on the surface of a tentatively designed structure are 
analyzed; a Process 2 in which an artificial high-cycle 
component is superimposed on the low-cycle compo- 
nent of the temperature fluctuations; numerical calcula- 



tions of accurate solutions for the motion equation and 
the energy equation for fluids, with the new temperature 
fluctuations obtained from Process 2 serving as bound- 
ary conditions; and a Process 3 in which the high-cycle 
components of temperature fluctuations occurring at 
structure surfaces are analyzed. Process 1 and Process 
3 are performed by a numerical calculation apparatus 
equipped with a control section to control time step siz- 
es. 
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Description 

BACKGROUND OF THE INVENTION 

FIELD OF THE INVENTION 5 

This invention relates to a thermal design method 
for structures and the optimum numerical calculation de- 
vices for such designs, and in particular to a design 
method and devices that take into account the effects 
of the temperature fluctuations of fluids that come into 
contact with the structures. 

RELATED ART 

Temperature fluctuations occur in areas where flu- 
ids of different temperatures are mixed together. If a 
structure exists in this area, repeated expansion and 
contraction of the structure surface due to the tempera- 
ture fluctuations will occur. In an airplane engine, for in- 
stance, the materials of the mechanisms which come 
into contact with heated gases and coolants are subject 
to these effects over long spans of time. 

Due to the expansion and contraction, irregularly 
fluctuating thermal stresses in the structure will occur. 
Since thermal stress usually occurs in the elastic parts 
of the structure, the frequency of expansion and con- 
traction can attain orders of 10 7 and above, resulting in 
so-called high-cycle thermal fatigue. Since cracks in the 
structure surface may develop depending on the ampli- 
tude, period, and frequency of the temperature fluctua- 
tions, high-cycle thermal fatigue has been considered 
to be an important design factor from the perspective of 
ensuring structural safety. 

The problem of high-cycle thermal fatigue is partic- 
ularly serious in fields where high levels of safety are 
required in a structure. For instance, in a fast reactor for 
a nuclear power plant, since high-temperature fuel 
channels and low-temperature control rod channels are 
mixed at relatively high velocities in the vicinity of the 
core exit, an important item in the course of develop- 
ment is structural design that takes into account the ef- 
fects of the turbulence that occurs in such areas. 

When considering the effects of heat in structural 
design, the velocity of heat flux is an important factor. 
Where the change in thermal flux over time is small, the 
orthodox engineering approach applies the time-aver- 
aged transport model with satisfactory results. The mod- 
el solves for time-averaged approximate solutions (i.e., 
the effective values for the state variables of a fluid) to 
the Navier-Stokes equation, the motion equation for flu- 
ids, and the energy equation for fluids, and describes 
only the low-cycle component of thermal flux behavior. 
Therefore, when the model is employed to analyze high- 
cycle components, engineering models such as the tur- 
bulence model are introduced, but since the solution is 
averaged overtime, it is not possible to analyze with pre- 
cision the high-cycle components with frequencies be- 



yond a certain range. Thus, in practice, the model is not 
applicable to the fast reactor design mentioned above. 

Further, solving for an accurate solution instead of 
such an approximate solution has also been consid- 
ered. However, the direct simulations necessary to ob- 
tain an accurate solution would require a tremendous 
amount of computing time since numerical calculations 
would be necessary which were based on a mesh ar- 
rangement fine enough to describe the minute struc- 
tures in the flow area and which employed time step siz- 
es small enough to describe the behavior of the random- 
ly fluctuating fluid. Therefore, it is difficult even with the 
aid of today's fastest computers to analyze fluid temper- 
ature fluctuations within the entire space of a fast reac- 
tor. A lurther, crucial problem is that even if such a 
lengthy computation was to be executed, as long as cur- 
rent methods for direct simulation are to be applied, 
there will be cases where the final solution cannot be 
attained due to the occurrence of divergent or unstable 
numerical values during calculations. This happens be- 
cause while the behavior of the state variables in the 
flow area should be described in terms of physical vi- 
brations, physically meaningless numerical errors occur 
in the computational process, since although direct sim- 
ulation is a method that solves for accurate solutions, 
mesh calculations performed in the process involve ap- 
proximating stages which are characteristic to comput- 
ers. 

In short, currently, with respect to structures which 
are subject to the effects of large fluid temperature fluc- 
tuations over time, analytical methods for highly precise 
structural design are practically non exist. As a result, 
when newly designing complex structures such as fast 
reactors, experiments using full size mockups are held 
to be necessary. As long as numerical analysis cannot 
be performed, there is no choice but to make evalua- 
tions on the basis of direct and reliable experiments. 

Problems related to the orthodox technology men- 
tioned above are summarized below from the following 
three perspectives: (1 ) deriving an approximate solution 
from the time-averaged transport model; (2) deriving an 
accurate solution from direct simulation; and (3) evalu- 
ation by means of mockup experiments. 

1. Deriving an approximate solution from the time- 
averaged transport model. 

Since this method merely establishes a turbulence 
model and solves for a time-averaged solution, it cannot 
describe the high-cycle component of fluid temperature 
fluctuations with precision. In areas of flow turbulence, 
it is mainly the high -cycle component of thermal fluctu- 
ations that affects high -cycle thermal fatigue. Therefore, 
the use of a model is meaningful only if the model can 
afford a high degree of precision. Further, even if it's use 
on high-cycle thermal fatigue analysis was restricted 
only to areas with slow thermal flux, as in the orthodox 
approach, current numerical calculation methods con- 
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sume too much time. Moreover, since the time step siz- 
es are inappropriate, it is not always possible to arrive 
at a solution, a problem also encountered in the direct 
simulation method to be discussed next. 

s 

2. Deriving an accurate solution from direct simulation. 

A tremendous amount of time is required. There- 
fore, this method cannot be applied to structures beyond 
a certain size and degree of complexity. Even if the prob- 10 
lem of computing time is tolerated, there will be cases 
in which solutions cannot be attained due to the diver- 
gence or instability of numerical values. 

3. Evaluation by means of mockup experiments. 

The method is direct and will yield certain results. 
However, when experimental results are deemed unsat- 
isfactory, design revisions, the construction of a new 
mockup, and the gathering of new sets of data are nec- 
essary. Design revisions are owed largely to trial and 
error, and the decision on a final structure often con- 
sumes much time and effort. To design a fast reactor, 
for instance, several years of experimenting and several 
hundred million yen's worth of funding would be neces- 
sary. Further, in designing a very large structure, it is 
often difficult to construct full size mockups. Errors 
which go hand in hand with experiments is another prob- 
lem, and, naturally, any compromises made in the con- 
struction of the mockup will also compromise the relia- 
bility of the experimental data. 

SUMMARY OF THE INVENTION 

In view of the problems discussed above, the pur- 
poses of this invention are classified below into several 
steps: 

1. In deriving an approximate solution from the 
above-mentioned time-averaged transport model, 
to provide a numerical calculation device that 
reduces computing time. 

2. In deriving an accurate solution from the above- 
mentioned direct simulation, to provide a numerical 
calculation device that makes the reduction of com- 
puting time and the stabilization of numerical values 
possible. 

3. In any given area of flow, to provide a method of 
analysis particularly with respect to the thermal 
design of structures in areas of flow turbulence. 
Here, a design method that greatly reduces design 
evaluation time by jointly employing the numerical 
calculation devices of both (1 ) and (2) above is dis- 
closed. 

Thus preferred embodiments of the invention, with 
respect to complex structures affected by fluid temper- 
ature fluctuations, establish and disclose an analytical 



method for thermal design, and provide optimum nu- 
merical calculation devices for the implementation of 
said method. The purpose of these embodiments is to 
make it possible to discontinue mockup experiments, re- 
duce designing time and cost, eliminate uncertain fac- 
tors in experiments, grasp numerically the amount of 
leeway offered by a design, and accumulate know-how 
for the sake of future designs. 

In order to achieve this object, a currently preferred 
embodiment provides a thermal design method for 
structures. This method comprises numerically calculat- 
ing time-averaged solutions of a motion equation and 
an energy equation for fluids to primarily analyze a low- 
cycle components of fluctuations in fluid temperature, 
said fluctuations occurring within space of a tentatively 
designed structure, superimposing an artificial high-cy- 
cle component on the low-cycle component of the fluc- 
tuations in fluid temperature obtained by the primary 
analysis to obtain new fluid temperature fluctuations, 
and numerically calculating accurate solutions of a mo- 
tion equation and an energy equation for fluids, using 
the new fluid temperature fluctuations obtained by the 
superimposition as boundary conditions, to secondarily 
analyze a high-cycle components of the fluctuations in 
fluid temperature occurring at any points of interest with- 
in the space of the structure, thereby generating infor- 
mation on temperature fluctuations at structure surfaces 
so as to design a structure taking into account effects of 
fluid temperature fluctuations. 

In accordance with this method, based on informa- 
tion on the form of the tentatively designed structure ob- 
tained by the primary analysis step, numerical calcula- 
tions are performed to derive time-averaged solutions 
of the motion equation and the energy equation for fluids 
that come into contact with the structure. As a result, the 
low-cycle components of the fluid temperature fluctua- 
tions that occur within the space of the structure are an- 
alyzed. Next, in the superimposition step, an artificial 
high-cycle component is superimposed on the low-cycle 
component of the fluid temperature fluctuation. This su- 
perimposition is based on empirically established rules 
on the behavior of fluid temperature fluctuations. In the 
secondary analysis step, the new fluid temperature fluc- 
tuations thus derived are employed as boundary condi- 
tions, and numerical calculations are performed to ob- 
tain accurate solutions to the motion equation and the 
energy equation for fluids, and the high-cycle compo- 
nent of the fluid temperature fluctuations occurring at 
each point of interest within the space of the structure 
is analyzed. Through these steps, information on the 
thermal fluctuations at structure surfaces necessary for 
the evaluation of the thermal fatigue occurring in struc- 
tures can be generated with precision. 

The above-mentioned superimposition step in- 
cludes superimposition of a high-cycle component arti- 
ficially generated from a series of random digits on the 
low-cycle component of the fluctuations in fluid temper- 
ature obtained by the primary analysis step. The 
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grounds for this is the fact that high-cycle components 
are white noise, and thus the superimposition makes it 
possible to compute the high -cycle components with 
precision. 

In another aspect of the invention, a numerical cal- $ 
culation apparatus for analyzing fluid temperature fluc- 
tuations that occur within space of structures. This ap- 
paratus comprises storage sections for storing a turbu- 
lence model approximately expressing a behavior of tur- 
bulence as well as storing a motion equation and an en- 
ergy equation for fluids which govern the behavior of tur- 
bulence, calculation sections for performing iterative 
computation of state variables of a flow area, in accord- 
ance with the turbulence model and the motion equation 
and the energy equation for fluids stored in the storage 
section, the iterative computation providing a new state 
variable based on the current state variable, with a given 
time step size between the new state and the current 
state, and a control section for controlling said time step 
size, in accordance with the state variables calculated 
by the calculation sections, thereby performing numer- 
ical calculations for time-averaged solutions to the mo- 
tion equation and the energy equation for fluids to ana- 
lyze the fluid temperature fluctuations. 

A turbulence model and the motion equation and 
the energy equation for fluids, the two equations gov- 
erning the behavior of fluids, are temporarily stored in 
the storage section. Next, in accordance with the turbu- 
lence model and the motion equation and the energy 
equation for fluids, iterative computation of the state var- 
iables of a flow area are performed by the computing 
section, the iterative computation providing a new state 
variable based on the current state variable, with a given 
time step size between the new state and the current 
state. Here, in accordance with the behavior of the state 
variable which is continuously changing, the time step 
sizes are controlled in order to optimize the computing 
time as well as the stability of numerical values during 
computation. 

In another aspect of the invention, a numerical cal- 
culation apparatus for analyzing fluid temperature fluc- 
tuations that occur within space of structures. This ap- 
paratus comprises storage sections for storing a motion 
equation and an energy equation for fluids which govern 
the behavior of turbulence; calculation sections for per- 
forming numerical calculations, through direct simula- 
tion, to solve the motion equation and the energy equa- 
tion which have been stored in the storage section, and 
based on calculation results, for performing iterative 
computation of the state variables of a flow area to pro- 
vide a new state variable based on the current state var- 
iable, with a given time step size between the new state 
and the current state, and a control section for control- 
ling the time step size, in accordance with the state var- 
iables calculated by the calculation sections, thereby 
performing numerical calculations for accurate solutions 
to the motion equation and the energy equation for fluids 
to analyze the fluid temperature fluctuations. 
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The motion equation and the energy equation for 
fluids which govern the behavior of flow areas are stored 
temporarily in the storage section. Subsequently, in ac- 
cordance with these equations, iterative computation of 
the state variables of a flow area are performed, through 
direct simulation, by the computing section, the iterative 
computation providing a new state variable based on the 
current state variable, with a given time step size be- 
tween the new state and the current state. Here, in ac- 
cordance with the behavior of the state variable which 
is continuously changing, the time step sizes are con- 
trolled in order to optimize the computing time as well 
as the stability of numerical values during computation. 

The above-mentioned control section may com- 
prise a fuzzy control section for controlling the time step 
size in accordance with the relative variation of the state 
variables calculated by the calculation section, an infer- 
ence section for inferring in advance locations of areas 
in which physically meaningful vibrations may occur, us- 
ing general rules related to fluid behavior, and an adjust- 
ment section for adjusting the time step size, only when 
vibrations appear in state variables during the calcula- 
tions for those areas which have not been inferred by 
the inference section as the areas in which the physi- 
cally meaningful vibrations may occur. As a result, the 
elimination of physically meaningless vibration compo- 
nents can be attempted. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 Block diagram showing the components of 
the numerical calculation device related to EMBODI- 
MENT!. 

Fig. 2 Diagram showing the membership function 
used to evaluate by means of a fuzzy set the size of the 
relative variation of state variable 

Fig. 3 Diagram showing the membership function 
used to evaluate the size of relaxation factor a. 

Fig. 4 Diagram showing the conversion of the size 
of relative variation to numerical stability. 

Fig. 5 Diagram showing the relationship between 
relative variation and a 

Fig. 6 Diagram showing the sliding of the member- 
ship function shown in Fig. 2. 

Fig. 7 Diagram showing the time taken to obtain a 
time-averaged solution for region in the vicinity of the 
fast reactor's core exit. 

Fig. 8 Block diagram of the numerical calculation 
device of EMBODIMENT 2. 

Fig. 9 Diagram showing the membership function 
for evaluating distance from a solid wall. 

Fig. 10 Diagram showing an area used in inference. 

Fig. 11 Diagram showing the effect of a solid wall 
expressed by converting the weighted distance from the 
solid wall. 

Fig. 12 Diagram showing temperature fluctuations 
occurring in sodium fuel in the mixing flow area. 

Fig. 1 3 Flowchart for the design method of EMBOD- 
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I ME NT 3. 

DESCRIPTION OF THE PREFERRED 
EMBODIMENTS 

EMBODIMENT 1 

The EMBODIMENT relates to a numerical calcula- 
tion device that shortens the computing time for the 
above-mentioned time -averaged transport model. 

Figure 1 shows a block diagram of the components 
of the numerical calculation device. The components of 
said device will be explained using this diagram. 

As shown in the diagram, said device has a momen- 
tum calculation section 2, a turbulence calculation sec- 
tion 4, and an energy calculation section 6 which respec- 
tively calculate a fluid's state variables, namely, momen- 
tum, turbulence and energy, at any given moment. What 
is notable here is that since said device approximates 
the effects of turbulence, the turbulence calculation sec- 
tion exists independent of the momentum calculation 
system. In other words, if an accurate solution instead 
of an approximate solution is required, the turbulence 
calculation will be inseparable from the momentum cal- 
culation section 2. 

The state variables calculated by the three calcula- 
tion sections are represented by the character <f>, and 
the state variable at time tn is written as 4>n. Thus, said 
device can be described as a device that takes an input 
<t>n of a given moment and calculates for <|>n+1, then fur- 
ther calculates for (j>n+2, in loop fashion. 

The state variable calculations performed by said 
device are in accordance with the Navier-Stokes equa- 
tion expanded for numerical analysis and the energy 
equation for fluids, etc. In order to derive a time-aver- 
aged solution that takes into account the effects of the 
high-cycle component, a turbulence model is used. Not 
only the k-e model, the linear time model generally used 
in orthodox engineering, but the algebraic stress model 
or the Reynolds stress model with their more precise 
quadratic approximations can also be used. 

The above describes the components of the calcu- 
lation section of said device. However, the characteris- 
tics of said embodiment is in the control section 8 which 
controls the value o1 the time step size At of the calcu- 
lation from moment to moment. 

The control section 8 keeps constant watch over 
A(|>/<{>, the relative variation of the state variables, and is 
a type of fuzzy controller which maintains the optimum 
time step size At. Here, A<|> is defined as: 

Att>n=(<|>n+1- (|>n)/<|>n, 
and At is determined using the equation 
At=a -Ate. 

Ate is a fixed value that serves as a basis for the 
time step size, and a is called a relaxation factor, which 
is a control parameter falling within the range from 0 to 
1 that is determined by the control section 8 . The control 



section B refers to the following numerical values in de- 
termining relaxation factor a: 

a. from the momentum calculation section 2, the rel- 
5 ative variation of the momentum in directions x, y, 

and z of the fluid, respectively: 

Au/u, Av/v, Aw/w; 

b. from the turbulence calculation section 4, the 
kinetic energy of the turbulence, the rate of dissipa- 

10 tion of kinetic energy, and the relative variation of 
the amplitude of the temperature fluctuations, 
respectively: 

Ak/k, Ae/e, AG* 2 /9" 2 
15 c. from the energy calculation section 6, the relative 
variation of enthalpy: 

Ah/h. 

When the relative variation A$/$ of all the above fac- 

20 tors converge on 0, each state variable 4> is solved. 

Although the main function of the control section 8 
is to control relaxation factor a, when the convergence 
of the relative variation or the numerical stability is poor 
in spite of the controlling of relaxation factor a, an adap- 

2S tive function is also provided in which the membership 
function itself which is used by the fuzzy controller is 
adjusted, as explained below. 

The functions of the above-mentioned construction 
are explained below. 

30 First, the state variable (fn at a given moment is in- 
putted to the momentum calculation section 2. The mo- 
mentum calculation section 2, the turbulence calculation 
section 4, and the energy calculation section 6 perform 
calculations to obtain a time-averaged approximate so- 

3S lution, outputting the state variable (|>n+-1 . <|>n+1 is again 
inputted to the calculation section, and the state variable 
calculations are performed moment by moment. In the 
process, the control section 8 evaluates the numerical 
stability of the state variables on the basis of the relative 

4Q variation of each state variable, and controls the relax- 
ation factor a on the basis of the following rules: 

(a) If the relative variation is negative and very large, 
change a to a very large value. 
4S (b) If the relative variation is negative and large, 
change a to a large value. 

(c) If the relative variation is negative and medium, 
maintain the present a value. 

(d) If the relative variation is positive and large, 
so change a to a small value. 

(e) If the relative variation is positive and very large, 
change a to a very small value. 

These rules were determined on the basis of the 
55 following empirically established rules. Relative varia- 
tion is large when the numerical value is far from the 
solution, and during this time a should be kept small so 
as to deal with the large variation. On the other hand, 
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as the relative variation approaches the solution, the 
convergence value can be reached quickly with large 
time step sizes. 

Though this set of rules can be implemented with 
ease when there is only one state variable, since the 
solution that is being sought involves several state var- 
iables, a certain number ol analytic processes are nec- 
essary to determine a from the several relative varia- 
tions. The processes are explained with the help of di- 
agrams. 

Figure 2 shows a membership function which is 
used to evaluate by means of fuzzy convergence the 
relative variation of a state variable <|>. Relative variation 
is plotted on the horizontal axis, and the grade is plotted 
on the vertical axis. The horizontal axis is divided into 
five sections in a Gaussian profile, with the sets for 
"small," "small medium," "medium," "medium large," 
and "large" plotted respectively from left to right. Figure 
3 shows the membership function used to evaluate the 
size of relaxation factor a, with a on the horizontal axis 
and the grade on the vertical axis. The horizontal axis 
is divided into five sections in a Gaussian profile, with 
the sets for "quite small," "somewhat small," "keep," 
"somewhat large," and "quite large" plotted respectively 
from left to right. 

For instance, if only two relative variations Au/u and 
Aw/w are considered, given the following conditions: 

Au/u: the relative variation is negative and medium 
(with a grade of 0.6) 

Aw/w: the relative variation is positive and quite 
large (with a grade of 0.8), 
then the former yields "keep the present a," while 
the latter yields "set a at a quite small value." 

Figure 4 shows a graph in which the size of the rel- 
ative variation has been converted to express numerical 
stability, and shows the loci of Au/u and Aw/w. If relative 
variation is medium, the current a should be kept, while 
if relative variation is large, then small a should be ap- 
plied. The relationship between the relative variations 
and a is shown in Figure 5. Here, it should be noted that 
in making the conversion from Figure 4 to Figure 5, the 
respective grades were weighted. In Figure 5, a is de- 
termined as follows: The center of gravity G of the area 
formed by the union of the two areas of the "keep" set 
and the "small" set in Figure 5 is calculated and project- 
ed onto the base set. The resulting point P gives the 
value for a. 

In this example, the explanation involves two rela- 
tive variations. Cases with three or more relative varia- 
tions can be solved in a similar manner. In any event, 
by controlling a with the control section 8, not only is 
considerable time reduction possible in comparison to 
prior methods, but problems such as numerical diver- 
gence encountered with prior methods that base their 
calculations on a fixed time step size Ate have been 
practically eliminated. 
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Next, a second function of the control section B will 
be explained. 

By optimum control of relaxation factor a, the con- 
trol section 8 aims to facilitate numerical analysis. How- 
5 ever, when numerical stability or convergence is unsat- 
isfactory despite the control, the adaptive function of the 
control section 8 is used. This function keeps a constant 
watch over the slope of the relative variation, and when 
the improvement of the slope falls below a certain 
10 threshold, the membership function shown in Figure 2 
is slid for the base set. Figure 6 shows the result of slid- 
ing the membership function of Figure 2. Here the con- 
trol parameter is the amount of sliding. For each amount 
slid, the control section 8 records the degree of improve- 
rs ment in the slope of the relative variation, and optimizes 
the slide amount. In other words, the control section 8 
is also equipped with an adaptive function. This con- 
cludes the summary of the numerical calculation device 
related to EMBODIMENT 1 . 
20 Several mathematical models can be used as the 
calculation method for the calculation section of said de- 
vice, such as: 

(a) numerical calculation method: finite difference 
25 method 

(b) time integral method: revised ICE method 

(c) convective term evaluation: QUICK method 

The actual computing time taken by a supercom- 
30 puter to find the time-averaged solution lor the area in 
the vicinity of the core exit of a fast reactor using these 
calculation methods is shown in Figure 7. In the graph, 
curve A shows the time needed to solve for different ct's 
(fixed value) using said device without resorting to the 
35 functions of the control section 8, whereas curve B 
shows the time needed when the functions of the control 
section 8 were applied. Curve B is shown as a straight 
line because a ranges from 0 to 1 . While from seven to 
sixteen hours were needed without the control section 
40 8, use of the control section 8 resulted in an improve- 
ment reducing the hours needed to approximately 3 
hours. 

EMBODIMENT 2 

45 

A numerical calculation device which, like EMBOD- 
IMENT 1 , makes it possible to reduce computing time 
and derive a stable solution is explained. The compo- 
nents of EMBODIMENT 2 are similar to that of EMBOD- 
so IMENT 1 ; however, EMBODIMENT 2 is a direct simula- 
tion device that solves not for a time-averaged solution 
but for an accurate solution, and therefore differences 
to that effect exist. 

Figure 8 is a block diagram that shows the compo- 
se nents of the numerical calculation device of EMBODI- 
MENT 2. Using this diagram, the differences in the com- 
ponents of said device and those of the device of EM- 
BODIMENT 1 will be explained in main. Components 
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common to both devices have been marked with the 
same symbol. 

As shown in the diagram, said device also has a 
momentum calculation section 2 and an energy calcu- 
lation section 6 to calculate the respective state varia- $ 
bles of the fluid, namely, the momentum and the energy, 
at a given moment, . However, the turbulence calcula- 
tion section 4 that was necessary in EMBODIMENT 1 
does not exist here. This is because since said device 
treats the effects of turbulence not by approximation, but 
by directly solving the Navier-Stokes equation for each 
time and location, the turbulence calculations cannot be 
separated from the momentum calculation section 2. 
Thus, a turbulence model is also unnecessary. 

Further, said device has a control section 10 with 
functions which differ from EMBODIMENT 1. The con- 
trol section 8 of EMBODIMENT 1 has as its functions 
the controlling of relative variation a and the sliding of 
membership functions; in contrast, the control section 
10 of said device uses empirically established rules and 
common sense knowledge as rules for controlling rela- 
tive variation a. 

The operations of the components described above 
are explained below. 

The basic operations of said device are the same 
as those for the device in EMBODIMENT 1. However, 
this device deals with the physical phenomena of ran- 
dom high-cycle vibrations, and must therefore be con- 
sidered differently from EMBODIMENT 1. In other 
words, although said device solves for an accurate so- 
lution, as long as the calculations are made for mesh 
cells with finite sizes, approximations or errors will occur 
Such errors have undesirable effects that cannot be ig- 
nored in the analysis of turbulence. This is because mi- 
croscopic high-cycle vibrations in the turbulent area, 
that is physically meaningful vibrations (physical vibra- 
tions) are joined by physically meaningless vibrations 
caused by calculation errors (numerical vibrations), and 
in a worst case scenario, divergence of numerical val- 
ues may occur during calculations. 

Thus, the control section 10 of said device makes 
an initial qualitative inference to determine the possible 
flow geometries given the shape of a certain structure 
and predict the areas where physical vibrations might 
possibly occur. Vibrations occurring in areas other than 
the anticipated areas are judged to be numerical vibra- 
tions, and attempts are made to eliminate them by ad- 
justing relaxation factor a. 

This function of control section 10 is based on em- 
pirically established rules and common -sense knowl- 
edge in hydrodynamics. The four common-sense rules 
programmed into control section 1 0 are: 

1. Fluid flows from inlet to outlet. 

2. Fluid flows along a solid wall. 

3. The normal velocity to the solid wall is zero. 

4. Momentum changes its direction in the presence 
of obstacles in the downstream direction. 



Control section 1 0 uses the four rules to make qual- 
itative inferences about the fluid, and to guess the pos- 
sible areas where physical vibrations might occur. 

An instance of qualitative inference is explained us- 
ing a diagram. 

Figure 9 shows a membership function used to eval- 
uate the distances from solid walls. The horizontal axis 
shows the distances from the solid walls, and the vertical 
axis shows the grade. The sets "near," "medium, " and 
"tar'' are plotted on the horizontal axis according to their 
distances from the solid wall, starting with "near." Figure 
10 shows the area for the inference. Here the point Q 
which is at distances D1 from the left wall and D2 from 
the upper wall will be examined for the effects of the solid 
walls. 

Based on the above assumptions, points D1 and 
D2 are judged as follows from Figure 9: 

(a) D1: Medium, with grade of 0.1 

(b) D2: Near, with grade of 0.9 

Based on this, the distances from the solid walls are 
weighted and converted to express the degree of effect 
from the solid walls, giving the graph in Figure 11. In 
Figure 11, the sets "small," "medium," and "large" are 
plotted in order of size of effect, starting with "small." As 
in EMBODIMENT 1 , the center of gravity G for the shad- 
ed area in Figure 1 1 is calculated and projected onto the 
base set. The point P thus obtained gives the effect of 
the solid walls. 

Next, starting from the upstream side, the flow di- 
rection in each calculation cell is determined and weight- 
ed with the wall effect from Figure 11 , and the vorticity 
in each cell is obtained. Based on the vorticity, the areas 
where physical vibrations may occur are inferred. Vibra- 
tions occurring where physical vibrations cannot possi- 
bly occur are deemed to be numerical vibrations, and 
the relaxation factor a is changed in order to eliminate 
the error component. 

in this example : an area bounded by two solid walls 
is explained. For areas bounded by three or more solid 
walls, the same approach can be applied. In any case, 
through the qualitative inference of control section 10, 
the problem of numerical divergence experienced in pri- 
or devices is solved. 

This concludes the summary of the numerical cal- 
culation device related to EMBODIMENT 2. 

Several mathematical models can be used as the 
calculation method for the calculation section of said de- 
vice, such as: 

(a) The numerical calculation method: the finite dif- 
ference method 

(b) The time integral method: the Leap-Frog method 

(c) Convection term evaluation: the tertiary wind- 
ward finite difference method 

By employing these calculation methods in said de- 
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vice, engineering problems involving the behavior of 
high-cycle temperature fluctuations, which were thought 
to be impossible to solve, can be evaluated with a high 
degree of accuracy using existent computing resources. 

EMBODIMENT 3 

EMBODIMENT 3 of the invention is explained be- 
low. 

EMBODIMENT 3 makes use of the devices de- 
scribed in EMBODIMENT 1 and EMBODIMENT 2 and 
relates to a method for the analytic design of structures 
in any given flow area. First, in order to understand the 
process flow of EMBODIMENT 3, the empirical facts re- 
lated to temperature fluctuations are explained. 

In general, when fluids of different temperatures are 
mixed, the fluid temperature fluctuations in the mixing 
area have the following characteristics: 

Characteristic 1 : The low-cycle component of the 
fluid temperature fluctuations is dominant. 
Characteristic 2: The statistical characteristics of the 
high-cycle component of the fluid temperature fluctua- 
tions conform to those for white noise. 

For instance, in the vicinity of the core exit of a fast 
reactor, high -temperature sodium fuel and low-temper- 
ature coolant mix together. The flow velocity of each fluid 
is approximately 7 meters per second. In this case, past 
experimental data indicates that the contribution of vi- 
brations in the order of 0 to 20 hertz to the fluid temper- 
ature fluctuations occurring in the mixing flow area is 
extremely large, whereas the contribution of vibrations 
exceeding 20 hertz is extremely small and remains at a 
constant level. 

Temperature fluctuations occurring in the sodium 
fuel at the mixing flow area is explained using Figure 1 2. 
In the figure, the horizontal axis shows time, and the ver- 
tical axis shows temperature. Curve A is a description 
of actual temperature observations, and curve B was 
obtained by averaging the component in A which was 
causing the fine vibrations, namely, the high -cycle com- 
ponent, and bringing out the low-cycle component which 
forms the basis of temperature fluctuations. X is a spec- 
imen of the high -cycle component's wavelength, and the 
frequency can be determined from its reciprocal. W is 
an example of the high-cycle component's amplitude. 

As shown in the figure, although the actual fluid tem- 
perature fluctuations appear to be completely irregular 
in behavior, on the whole, the low-cycle component 
shown by curve B is dominant, and the shape of the fluc- 
tuations is the result of superimposing the random high- 
cycle component, or white noise, on the low-cycle com- 
ponent. This high-cycle component occurs completely 
randomly, reaching an order of several kilohertz. This 
illustrates the above-mentioned characteristics 1 and 2. 

From characteristics 1 and 2, it can be seen that 
fluid temperature fluctuations can be analyzed by per- 
forming calculations separately for the low-cycle com- 
ponents and the high -cycle components of the fluctua- 



tions. Therefore, in EMBODIMENT 3, the low-cycle 
component is described as a time-averaged solution, to 
which a random high-cycle component is artificially su- 
perimposed. The synthesized vibrations are used as 
boundary conditbns for recalculation, with an aim to ob- 
tain an accurate solutbn. Here, the optimum methods 
for obtaining a time-averaged solution and an accurate 
solution are provided by the numerical calculation de- 
vices of EMBODIMENT 1 and EMBODIMENT 2, re- 
spectively. 

Based on the above discussion, each process in the 
design method of EMBODIMENT 3 is explained. The 
structural design of a fast reactor will be used as an ex- 
ample to explain said embodiment. Figure 13 is a flow- 
chart of the design method. The design method consists 
of the following procedures. 

Procedure 1: Determining the characteristics of the 
low-cycle temperature fluctuations. 

First, a tentative design for the structure is made. 
Next, the time-averaged solution for the Navier-Stokes 
equation is calculated using EMBODIMENT 1 or an 
equivalent device, yielding the low-cycle component of 
the fluid temperature fluctuations occurring within the 
space of the tentatively designed structure. 

Procedure 2: Super imposition of the high -cycle 
component. An artificial high -cycle component is super- 
imposed on the low-cycle component of the fluid tem- 
perature fluctuations obtained in Procedure 1 . The ran- 
domness of the high-cycle component is imitated using 
arithmetical random digits. In other words, both the am- 
plitude and frequency of the high-cycle component are 
generated using random numbers. The upper limit for 
the frequency is adjusted to several kilohertz. The result 
is a new model that can be used in Procedure 3. 

Procedure 3: Determining the characteristics of the 
fluid temperature fluctuations. 

By means of EMBODIMENT 2 or an equivalent de- 
vice, numerical calculations are performed in order to 
obtain an accurate solution to the Navier-Stokes equa- 
tion, using the fluid temperature fluctuations of the mod- 
el obtained from Procedure 2 as boundary conditions. 
As a result, fluid temperature fluctuations that accurately 
reflect the effects of the high-cycle component can be 
obtained. 

The above is the design method related to EMBOD- 
IMENT 3. When actually designing a structure, the re- 
sults obtained from Procedure 3 are used as a basis for 
the calculation of temperature fluctuations on the struc- 
ture surface, and the thermal stress occurring as a result 
of the temperature fluctuations is used to evaluate ther- 
mal fatigue. When the results of the evaluation are not 
desirable, the shape of the structure is modified, and the 
process proceeds along route a in the figure to repeat 
procedures 1 through 3, thus evaluating the suitability 
of the structural design. However, in practice, by making 
a preliminary evaluation after the completion of Proce- 
dure 1 (route b in the figure), design time can be re- 
duced. In other words, based on the approximate solu- 
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tion of Procedure 1, locations on the structure surface 
where the occurrence of thermal fatigue is predicted can 
be selected, and procedures 2 and 3 can be executed 
for these selected spots only. This preliminary evalua- 
tion is appropriate to a certain extent given the fact that s 
low-cycle components are dominant in fluid temperature 
fluctuations and the fact that high-cycle components oc- 
cur randomly. In other words, the invention makes time 
reduction possible by analyzing empirical facts about 
fluid temperature fluctuations and developing a method to 
that distinguishes between the tow-cycle component 
and the high-cycle component 

EMBODIMENT 3 deals with both the low-cycle and 
high-cycle components. The reason why this is indis- 
pensable for high-precision analysis is explained. is 
Although the low-cycle component dominates in fluid 
temperature fluctuations, the effects of the high-cycle 
component are greater in high-cycle thermal fatigue. 
Further, the transfer of energy from the low-cycle com- 
ponent to the high-cycle component must also be con- 20 
sidered. Therefore, an analysis involving only one of the 
two components would be incomplete. Prior design 
methods only solved for a time-averaged solution based 
on a turbulence model, and were thus unable to gener- 
ate accurate information necessary for the evaluation of 2s 
high-cycle thermal fatigue. The accurate information 
necessary is provided for the first time by the design 
method of EMBODIMENT 3. 

The explanation of EMBODIMENT 3 refers to the 
design of a fast reactor as an example; however, the 30 
design method naturally is not restricted to fast reactor 
design. The invention has an extremely wide range of 
applications, ranging from large-scale structures includ- 
ing thermoelectric power plants, chemical plants, steel 
manufacturing plants, space structures, and various 35 
other thermal facilities to electronic elements that are 
affected by temperature fluctuations. 

As has been explained in detail above, the thermal 
design method for structures related to the invention, by 
taking note of the fact that the low-cycle components 40 
and the high-cycle components of fluid temperature fluc- 
tuations can be separated, makes it possible to analyze 
the high-cycle components of fluid temperature fluctua- 
tions, which could not be accurately calculated with prior 
methods. The procedures are divided in accordance *s 
with the separation of fluid temperature fluctuation com- 
ponents, first obtaining the low-cycle component, then 
superimposing an artificial high-cycle component in or- 
der to obtain an accurate solution, thus providing a 
method for the analysis of high-cycle components of so 
flow areas, where such methods can be said to have 
been previously non-existent in the field of engineering. 
Further, since a preliminary evaluation is possible after 
the low-cycle component has been obtained, the inven- 
tion paves the way to practical analysis by accurate so- ss 
lutions, a method that prior to this suffered serious draw- 
backs due to computing time. 

Further, the thermal design method for structures 



related to the invention focuses on the fact that the 
above-mentioned high-cycle component is white noise, 
disclosing an epoch-making method in which a series 
of random numbers are used to describe said compo- 
nent. Since this method of description makes it possible 
to treat the high-cycle component in an analytical man- 
ner, it allows for the realization of extremely superior de- 
sign reproducibility and reliability. 

Further, by using the two types of numerical calcu- 
lation devices related to the invention, owing to the op- 
timum control of time step sizes for the calculations, re- 
ductions in computing time, improvements in computing 
precision, and numerical stability can be realized in cal- 
culations for approximate and accurate solutions to the 
motion equation and the energy equation for fluids. 
Therefore, by introducing computing science into the 
field of thermal design where the practical use was 
thought to be difficult, the invention makes the design of 
complex and large-scale structures within the scope of 
existent computing capabilities possible. 

Thus, by a combination of the thermal design meth- 
od and the numerical calculation devices related to the 
invention, mockup experiments which were previously 
held to be necessary in structural design have been ren- 
dered unnecessary, and the optimization of designs can 
be realized with ease. For instance, the design of a fast 
reactor for a nuclear plant which usually requires several 
years of experimentation and several hundred million 
yen in funding for such experimentation can be evalu- 
ated within a time span of several months and at a cost 
of several million yen. Further, by eliminating the need 
for experimentation, and therefore eliminating the need 
to consider the inevitable elements of uncertainty and 
measurement errors, the reliability of designs can be im- 
proved. The invention has a wide range of engineering 
applications, and its effects are considerable. 

While there has been described what are at present 
considered to be preferred embodiments of the inven- 
tion, it will be understood that various modifications may 
be made thereto, and it is intended that the appended 
claims cover all such modifications as fall within the true 
scope of the invention. 



Claims 

1 . A thermal design method for structures comprising 
steps of: 

numerically calculating time-averaged solu- 
tions of a motion equation and an energy equa- 
tion lor fluids to primarily analyze a low-cycle 
component of fluctuations in fluid temperature, 
said fluctuations occurring within space of a 
tentatively designed structure; 
superimposing an artificial high-cycle compo- 
nent on the low-cycle component of the fluctu- 
ations in fluid temperature obtained by the pri- 
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increase the time step size when the relative varia- 
tion is smalt. 

6. The numerical calculation apparatus of claim 4, 
5 wherein the control section (8) appropriately cor- 
rects membership functions referenced by the fuzzy 
control section when degree of improvement in the 
slope of the relative variation is small. 

to 7. The numerical calculation apparatus of claim 4, 
wherein, when a plurality of state variables are to 
be calculated by the calculation section, the fuzzy 
control section weights the plurality of relative vari- 
ations related to the state variables and integrates 

15 the weighted relative variations, thereby controlling 
the time step size on the basis of the results of inte- 
gration. 

8. The numerical calculation apparatus of claim 7, 
20 wherein said fuzzy control section controls the time 
step size on the basis of the center of gravity of each 
element area of fuzzy sets, the area size of each 
element area having been corrected by multiplying 
the original element area, to which each of the rel- 
2S ative variations belongs, by a grade corresponding 
to each of the relative variations. 



mary analysis to obtain new fluid temperature 
fluctuations; and 

numerically calculating accurate solutions of a 
motion equation and an energy equation for flu- 
ids, using the new fluid temperature fluctua- 
tions obtained by the superimposition as 
boundary conditions, to secondarily analyze a 
high-cycle component of the fluctuations in fluid 
temperature occurring at any points of interest 
within the space of the structure, 
thereby generating information on temperature 
fluctuations at structure surfaces so as to 
design a structure taking into account effects of 
fluid temperature fluctuations. 

2. The thermal design method of claim 1 , wherein said 
superimposition step includes superimposition of a 
high-cycle component artificially generated from a 
series of random digits on the low-cycle component 
of the fluctuations in fluid temperature obtained by 
the primary analysis step. 

3. A numerical calculation apparatus for analyzing 
fluid temperature fluctuations that occur within 
space of structures, comprising: 

(a) storage sections (2) (4) (6) for storing a tur- 
bulence model approximately expressing a 
behavior of turbulence as well as storing a 
motion equation and an energy equation for flu- 
ids which govern the behavior of turbulence; 

(b) calculation sections (2) (4) (6) for perform- 
ing iterative computation of state variables of a 
flow area, in accordance with the turbulence 
model and the motion equation and the energy 
equation for fluids stored in the storage section, 
the iterative computation providing a new state 
variable based on the current state variable, 
with a given time step size between the new 
state and the current state; and 

(c) a control section (8) for controlling said time 
step size, in accordance with the state variables 
calculated by the calculation sections, 

thereby performing numerical calculations for time- 
averaged solutions to the motion equation and the 
energy equation for fluids to analyze the fluid tem- 
perature fluctuations. 

4. The numerical calculation apparatus of claim 3, 
wherein the control section (8) includes a fuzzy con- 
trol section which controls the time step size, in 
accordance with the relative variation of the state 
variable calculated by the calculation section. 

5. The numerical calculation apparatus of claim 4, 
wherein the fuzzy control section reduces the time 
step size when the relative variation is large, and 



9. The numerical calculation apparatus of claim 8, 
wherein the control section (8) modifies the mem- 

30 bership functions by sliding their positions on the 
base set. 

10. A numerical calculation apparatus for analyzing 
fluid temperature fluctuations that occur within 

35 space of structures, comprising: 

(a) storage sections (2) (6) for storing a motion 
equation and an energy equation for fluids 
which govern the behavior of turbulence; 
40 (b) calculation sections (2) (6) lor performing 

numerical calculations, through direct simula- 
tion, to solve the motion equation and the 
energy equation which have been stored in the 
storage section, and based on calculation 
45 results, for performing iterative computation of 

the state variables of a flow area to provide a 
new state variable based on the current state 
variable, with a given time step size between 
the new state and the current state; and 
50 (c) a control section (1 0) for controlling the time 

step size, in accordance with the state variables 
calculated by the calculation sections, 
thereby performing numerical calculations for 
accurate solutions to the motion equation and 
55 the energy equation for fluids to analyze the 

fluid temperature fluctuations. 

11. The numerical calculation apparatus of claim 10, 
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wherein said control section (10) comprises: 

a fuzzy control section tor controlling the time 
step size in accordance with the relative varia- 
tion ot the state variables calculated by the cal- s 
culation section; 

an inference section for inferring in advance 
locations of areas in which physically meaning- 
ful vibrations may occur, using general rules 
related to fluid behavior; and io 
an adjustment section for adjusting the time 
step size, only when vibrations appear in state 
variables during the calculations for those 
areas which have not been inferred by the infer- 
ence section as the areas in which the physi- is 
cally meaningful vibrations may occur. 

12. The numerical calculation apparatus of claim 11, 
wherein the inference section infers the locations of 
areas where the vibrations may occur, on the basis zo 
of the distance from a solid wall constituting a struc- 
ture. 

13. The numerical calculation apparatus of claim 11, 
wherein, in the presence of a plurality of walls con- 25 
stituting a structure, the inference section weights 

the respective distances from the walls and inte- 
grates the weighted distances to infer the locations 
of the areas where the vibrations may occur. 

30 

14. The numerical calculation apparatus of claim 13, 
wherein the inference section infers the locations of 
areas where vibrations may occur on the basis of 
the center of gravity of each element area of fuzzy 
sets, the area size of each element area having 35 
been corrected by multiplying the original element 
area, to which the respective distances belong, by 

a grade corresponding to each of the relative dis- 
tances. 



BNSDOCID: <EP 0709789A2_I_> 



11 



EP 0 709 789 A2 



8 



a 



0n 



CONTROL 
SECTION 



At=aAtc 



Au/u,Av/v,Aw/w 



1 



MOMENTUM 
CALCULATION 
SECTION 



TURBULENCE 
CALCULATION 
SECTION 



Ak/k 

Ae/e 

A6 ,2 /e'* 6 



ENERGY 
CALCULATION 
SECTION 



Ah/h 



<j>n+1 



Fig. 1 



SMALL SMALL MEDIUM MEDIUM MEDIUM LARGE LARGE 




RELATIVE VARIATION 



Fig. 2 



BNSDOCID: <EP 0709789A2J_> 



12 



EP 0 709 789 A2 



QUITE SOMEWHAT SOMEWHAT QUITE 

SMALL SMALL KEEP LARGE LARGE 




0.0 0.5 1.0 

RELAXATION FACTOR 



Fig. 3 



SMALL MEDIUM LARGE 




RELATIVE VARIATION 

Fig. 4 



13 

BNSDOCID: <EP 07O97B9A2J_> 



EP 0 709 789 A2 




Fig. 6 



14 



EP 0 709 789 A2 



TIME 




0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 
RELAXATION FACTOR a 



Fig. 7 



<!>n 



a 



10 



CONTROL 
SECTION 



At=aAtC Au/u,Av/v.Aw/w 





MOMENTUM 
CALCULATION 
SECTION 



ENERGY 
CALCULATION 
SECTION 



Ah/h 



Fig. 8 



BNSDOCID: <EP 07O97B9A2J_> 



15 



EP 0 709 789 A2 



NEAR MEDIUM FAR 




SOLID WALL 



Fig. 10 



16 



EP 0 709 789 A2 




B 



480.0 
475.0 
470.0 
465.0 
460.0 
455.0 
450.0 
*~ 445.0 
440.0 



o 

o 

LU 
DC 

< 

UJ 
Q. 



































1 w 






























J- 










i X 



30.0 32.0 



34.0 36.0 
TIME (sec.) 



38.0 40.0 



Fig. 12 



EP 0 709 789 A2 



LU 

CC 



CO 

■■■■ 

LL 



18 



BNSOOCID: <EP O7097B9A2J_> 



(19) 



Europaisches Patentamt 
European Patent Office 
Office europeen des brevets 



(12) 



(88) Date of publication A3: 

29.01.1997 Bulletin 1997/05 

(43) Date of publication A2: 

01.05.1996 Bulletin 1996/18 

(21) Application number: 95307699.9 

(22) Date of filing 27.10.1995 



( ii) EP 0 709 789 A3 

EUROPEAN PATENT APPLICATION 

(51) mtci A G06F 17/50, G21C 21/00 



(84) 


Designated Contracting States: 


(72) 


Inventor: Muramatsu, Toshiharu 




DE FR GB 




Higashiibarakl-gun, Ibaraki (JP) 


(30) 


Priority: 31.10.1994 JP 266773/94 


(74) 


Representative: 








Cross, Rupert Edward Blount et a I 


(71) 


Applicant: DORYOKURO KAKUNENRYO 




BO LILT WADE TENNANT 




KAIHATSU JIGYODAN 




27 Furnival Street 




Tokyo-to (JP) 




London EC4A1PQ (GB) 



CO 

< 

o> 

CO 

a> 
o 

o 

CL 
LU 



(54) A thermal design method and apparatus 

(57) A design method for structures is invented, tak- 
ing into ac count the effects of fluid temperature fluctu- 
ations, comprising: numerical calculations of time-aver- 
aged solutions for the motion equation and energy 
equation for fluids; a Process 1 in which the low-cycle 
components of the temperature fluctuations occurring 
on the surface of a tentatively designed structure are 
analyzed; a Process 2 in which an artificial high-cycle 
component is superimposed on the low-cycle compo- 



nent of the temperature fluctuations; numerical calcula- 
tions of accurate solutions for the motion equation and 
the energy equation for fluids, with the new temperature 
fluctuations obtained from Process 2 serving as bound- 
ary conditions; and a Process 3 in which the high-cycle 
components of temperature fluctuations occurring at 
structure surfaces are analyzed. Process 1 and Process 
3 are performed by a numerical calculation apparatus 
equipped with a control section to control time step siz- 
es. 
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